Macroscopic self trapping in BECs: analysis of a dynamical quantum phase transition 



O 
(N 



(N 



B. Julia-Di'az/ D. Dagnino/ M. Lewenstein,^' J. Martorell/ and A. PoUs^ 

^ Departament d'Estructura t Constituents de la Materia, 
Universitat de Barcelona, 08028 Barcelona, Spain 
^ ICREA-Institucio Catalana de Recerca i Estudis Avangats, Lluis Companys 23, 08010 Barcelona, Spam 
ICFO-Institut de Ciencies Fotdniques, 08860 Castelldefels (Barcelona), Spain 

(Dated: May 13, 2010) 

We consider a Bose-Einstein condensate in a double-well potential undergoing a dynamical tran- 
sition from the regime of Josephson oscillations to the regime of self-trapping. We analyze the 
statistical properties of the ground state (or the highest excited state) of the Hamiltonian in these 
two regimes for attractive (repulsive) interactions. We demonstrate that it is impossible to de- 
scribe the transition within the mean- field theory. In contrast, the transition proceeds through a 
strongly correlated delocalized state, with large quantum fluctuations, and spontaneous breaking of 
the symmetry. 
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I. INTRODUCTION 

Cold atom physics, e.g. Bose-Einstein condensates 
(BEC), provides a useful set up to tackle fundamental 
problems of quantum physics at a macroscopic scale. 
The physics of weakly interacting Bose gases has been 
described very successfully by employing a mean field 
approach, that is the Gross-Pitaevskii equation (GPE). 
Such approach implies that to a certain extent quantum 
fiuctuations and correlations are relatively unimportant 
in a large variety of phenomena, see for instance 



A. Quantum correlations and symmetry breaking 

The mean field approach can be improved by including 
in the theoretical description small (Gaussian) fiuctua- 
tions around the mean field solutions. These fiuctuations 
are described then by the Bogoliubov-de Gennes (BdG) 
equations, and have been intensively studied in experi- 
ments and theory Q . The approach based on BdG equa- 
tions has its limits also. In particular, it is completely in- 
adequate in strongly correlated quantum systems, such 
as for instance those that can be realized with ultracold 
gases in optical lattices, or in low dimensional gases (for 
recent reviews see [1, Interestingly, it has recently 
been pointed out by several authors that the BdG ap- 
proach breaks down in systems that undergo transition 
between two states that are well described by the mean 
field theory, but that differ dramatically in form, for in- 
stance having different symmetry properties. Transition 
between such states on the level of mean field descrip- 
tion is typically associated with dynamical instabilities 
of the GPE, and exponential growth and squeezing of 
the BdG fluctuation. The state of the system in such sit- 
uations becomes strongly correlated and entangled, and 
its description requires necessarily to go beyond the mean 
field theory. Several authors have pointed out this fact 
in various contexts: sonic analogues of "black holes" @, 
vortex nucleation in small atomic clouds [a-H, BECs in 



rotating ring superlattices [13], or in the ground state of 
BEC in a double well potential [T]| . 

Let us explain how the quantum correlations attain 
significance for the case of vortex nucleation in a rotat- 
ing BEC In this paper exact diagonalization meth- 
ods for small atomic clouds are used and compared with 
the mean field description, indicating directly the impos- 
sibility of a correct description of the transition within 
the mean field formalism. By performing the exact di- 
agonalization of the Hamiltonian of the system, the au- 
thors show that in the transition region the atoms are 
far from being in the same single particle state. Instead, 
two macro-occupied single particle states (MSPS) appear 
with comparable occupations. 

In this article we will show how the physics of BECs in 
a double-well resembles the above scenario. In the case 
of the rotating BEC, the control parameter driving the 
transition is the rotation frequency. In the case of the 
BEC in a double-well, it will be the interaction between 
the atoms, or, correspondingly, the height of the poten- 
tial barrier splitting the BEC. Contrary to the previous 
works, that focus on the properties of the ground state, 
or low energy states of the considered system, we deal 
here with the dynamical transition both for attractive 
and repulsive interactions, for which also the properties 
of the highest energy, or high energy states are relevant. 



B. BECs in double well potentials 

The physics of BECs on a double- well has attracted a 
great deal of attention since the first theoretical studies 
predicting the presence not only of Josephson-like oscil- 
lations of population, but also of self-trapped states, see 
Refs. Its recent experimental characterization by 

the Heidelberg group [13) opens the possibility of practi- 
cal applications of these bosonic junctions. 

As pointed out in Ref. [l|, most of the physics 
of bosonic Josephson junctions, including self-trapping, 
TT— modes and Josephson oscillations, can be understood 
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semiclassically. We will describe, however, how the 
transition from the Josephson-regime to the self-trapped 
regime is indeed dominated by quantum correlations. To 
that extent we will consider the exact diagonalization of 
the Bose-Hubbard (BH) Hamiltonian, and compare the 
results always to predictions of the semiclassical approach 
(in which quantum mechanical operators are replaced by 
their c- number averages). 

Our aim is to isolate the purely quantum aspects of 
the transition, defined as those which are either not un- 
derstandable in terms of a mean field description of the 
problem or, alternatively, by an elementary semiclassical 
analysis. 



C. The plan of the paper 

We will start, sections II-IV, discussing the properties 
of the BH Hamiltonian emphasizing the relevance of the 
condensed fractions, eigenvalues of the one-body density 
matrix, and the symmetry properties of the eigenstates 
of the Hamiltonian. By analyzing several quantities si- 
multaneously we are able to classify the different regimes 
into which the transition is divided and which are clearly 
beyond the mean field approximation taken previously by 
some authors. The quantum nature of the transition is 
clearly visualized. Then, in section V, we present the dy- 
namical consequences of the statical properties described 
before. The final section, VI, gives the conclusions of our 
analysis. 



II. DESCRIPTION OF THE MODEL 

As described in [l|, [l^l the two-mode Bose-Hubbard 
Hamiltonian for N atoms in the simple case of a double- 
well may be written as, 

H = {fiLifiL - ^) + nnifiR - I)) 

- J (a'j^aL + a^^aji^ - e{nL - ur) (1) 

where a\^\nL,nR) = ^/r^^~^T|nL + l,nR), aL\nL,nR) = 
y/rTElriL - l,nR), and [a|,aj] = Sij , i,j = L,R. L{R) 
stands for the left (right) side of the well. Our sign con- 
vention is such that U > 0{U < 0) corresponds to attrac- 
tive (repulsive) atom-atom interactions. We discuss both 
cases simultaneously as the technicalities are similar for 
both cases. We take J > 0. 

A small bias, < e <C J, is introduced which will 
ensure breaking of left-right symmetry. Numerically we 
consider, e/J— 10^^" in all the results presented, which 
intents to provide a realistic implementation of the sym- 
metric double well. 

The Hamiltonian is diagonalized in the + 1 di- 
mensional space spanned by the basis: {|A^, 0),|A^ — 
1,1),..., \1,N - l),|0,Ar)}. With, \nL,nR) = 



l/VnL!rrR!(a}^)"^(a]j)"^|vac). The dynamics of the sys- 
tem is governed by the ratio, A = NU/J, which is con- 
trolled either by varying the atom- atom interactions U, 
or, by changing the barrier height, J. 

Any single particle state may be written as, 
1*1 (61, 0)) = cos(6'/2)|L) -I- e''^sin(6i/2)|i?), a representa- 
tion usually employed to describe single qubits [l^. We 
define, \L{R)) = a^^^^jvac), and |± >= (1/V2)(|L) ± 
\R)). 

The most general A^-body state can be written as, 
I*) = Efe=i,Ar+ic(A)l^ + l-k,k-l). The number of 
atoms in each well for a given state is, A^^ = (^'jajja/jl^'), 
with P = L,R. The population imbalance of a state 
^' is defined as, z — {Nl — Nr)/N. The time evolu- 
tion of a state ^'(0) is governed by the time-dependent 
Schrodinger equation. 

To characterize the degree of condensation of the sys- 
tem at any given time we will make use of the one- 
body density matrix, p. For a state, *, we have, pij — 

(^'Ipyl^l') , with Pij = aluj , and i,j = L, R. The trace of 
p is normalized to the total number of atoms, TV. We will 
denote its two eigenvectors as i/'i(2) and its normalized 
eigenvalues (divided by the total number of atoms A) as 
ni(2), with ni > n2 > 0. We always have, ni+rt2 = 1. rii 
will correspond to the condensed fraction in the macro- 
occupied state ^pi ^. 

A measure of the spread of the state in the 
Fock basis will be given by the function 5* = 
~ Sfc^/ |c(^)P log(|c(A:)p). S is positive definite and 
has a maximum value for the equally populated state, 
c^(fc) = 1/(A^ -I- 1), and is zero for a maximally localized 
state, c{k) = Sk,ko- 

To clearly identify the region where the quantum fluc- 
tuations will play a major role we will always compare 
our numerical results to semiclassical predictions. The 
semiclassics of the problem may be obtained by replac- 
ing the creation and anhilation operators by c-numbers, 
dLiR) = ^'^L(_R)e*'''^<'^'. We can define the phase differ- 
ence as, ip = ifR — (fL- The semiclassical Hamiltonian 
reads, 

— = - Vl -z^cosip- —(Nz^ + N-2), (2) 

and the time evolution of the population imbalance and 
phase difference is obtained through the Heisenberg evo- 
lution equations, icll = — [i?, a^] and iclr = — aj^]. 
The resulting equations are, ^j- — — Vl — z"^ sin Lp, and 
— — 4^ + , ^ cos If , which were originally ob- 
tained by Smerzi et al. [ij] as the two-mode approxi- 



Other authors, e.g. Ref. [16L Il7l|. introduce the operators, Jx = 
{l/2)(a^a^ -I- aJjOi), Jy = (l/2i)[a\an - a^aj,), and = 
(1/2) (a^a^ —a'^an), which are related to p though, pi^^ = N/2+ 
{Jz), PRR = N/2-{J;,), plr = {Jx+iJy), and PR^ = (Jx-tJy)- 
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FIG. 1: Density plot of \cx{k)\^ , the horizontal (vertical) 
axes run through k — 1, . . . , A'^ + 1 (A = 1, . . . , A'^ + 1), for, 
starting from above, (left) NU/J = —4, —8, —12 and (right) 
NU/J = 0, 4, 8, 12. The left (right) column of each two 
displays the eigenvectors with positive (negative) imbalance. 
Blue corresponds to zero and white to the maximum value. 
N is 50. 
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Spectral index (i)/(N+l) 

FIG. 2: Spectral decomposition, \c\{k)\^ of the ground 
state (solid, black) and the highest state of the Hamilto- 
nian (dashed, red) in the basis spanned by (n_t,n_R) — 
{{N, 0), {N -1,1),..., (0, N)} for iV = 50. 11 U ^ -U the 
same plots would be obtained, but with the roles of the ground 
and highest excited states exchanged. 



mation to the Gross-Pitaevskii mean field description of 
the problem. 

Ref. defined the transition from tlie Josephson 
to the self-trapped regime as a dynamical feature, that 
is, a state initially prepared with a certain (z(0), (/?(0)) 
would undergo Josephson oscillations or would remain 
self-trapped on one side of the double-well if A is larger 
than a critical value. We will link that definition, which 
is a semiclassical one in our model, with the static prop- 
erties of the BH Hamiltonian. 



III. GENERAL PROPERTIES OF THE 
SPECTRUM 

Following the discussions of Refs. [l^, let us first 
present the spectral decomposition, |cA(fc)P, of the eigen- 
vectors of the Hamiltonian, H'i'x = E\'i/\, as a function 
of A, see Fig. [TJ Fig. [5] concentrates on the properties 
of the ground state and the highest excited state of the 
spectrum as we vary the value of A = NU /J. The main 
features relevant for our discussion are; 



A. Attractive interactions: ground state analysis 

When e = the ground state is symmetric 
and binomial for the non-interacting case, "^gs = 
(l/%/]V!)[(l/V2)(a[ -hajj)]^|vac). Then, as the interac- 
tion is increased the ground state becomes wider. Then, 
it becomes degenerate with the next excited state, and 



its distribution has eventually two differentiated peaks, 
thus becoming cat-like 18]. When e > 0, but small. 
Figs. [U and [2] show that, as we increase the interaction 
further, each of the formerly degenerate pairs develops 
a certain population imbalance and the distribution be- 
comes peaked on only one region of the Fock space. It 
is important to note that the left-right symmetry is bro- 
ken for a certain value of A, thus, the ground state will 
spontaneously acquire a large imbalance. Our figure dif- 
fers from Fig. 2 of Ref [l^, because we have artificially 
broken the symmetry with the help of a small bias. In 
this way the two degenerate eigenstates are distinguished 
according to their imbalance. Our figure emphasizes the 
symmetry breaking pattern, which occurs when the sys- 
tem evolves from cat-like to self-trapped. 



B. Repulsive interactions: analysis of the highest 
energy state 

The case of repulsive interactions is similar but the 
role played by the ground state is now played by the 
highest excited state of the Hamiltonian. With no inter- 
action, the highest eigenvector of the Hamiltonian is also 
binomial, ^-he = l/\/W.[{l/V2){al - a+)]^|vac), as the 
interaction increases, undergoes similar features as the 
ground state did in the attractive case. The ground state 
in this case, on the contrary, remains mostly binomial 
with a certain squeezing, as also discussed in Ref. [Toj . 
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FIG. 3: (a) Solid line: population imbalance z of the ground 
state. The shaded region corresponds to its dispersion, a^- 
The inset provides a logarithmic view of Oz, computed from 
the BH (dotted) and its semiclassical approximation (dashed) . 
(b) A measure of the spread of the ground state, S, nor- 
malized to its maximum value, which would correspond to 
a equidistributed state, (c) Solid and dotted lines depict 
the condensed fractions, n\ and ni, of the one body density 
matrix of the ground state of the Hamiltonian, with imbal- 
ance > as a function of A = NU/J. (d) The four ma- 
trix elements of p are depicted as a function of NU/J. The 
shaded region corresponds to the variance of the diagonal el- 
ements, aLL, where aij = yJ\{Pij) — {pij)^\ {o-rr turns out 

to be equal to aLL). The variance of the off-diagonal ones is 
negligible in this plot. The dashed lines in all the plots cor- 
respond to the semiclassical predictions (at e = 0): Zs.c. ~ 
if A < 2, Zs.c. = VI - 4/A2 if A > 2. pI';l = N/2 + N zs.o./2, 
p%y = N/2 - Nzs.,./2, pl% = p%% = NiVl-z'i.c.)/2. The 
number of atoms is, — 50. 



IV. ATTRACTIVE INTERACTIONS: 
QUANTUM FLUCTUATIONS 

To characterize the problem niore deeply, we now turn 
to other static properties of the eigenstates. In Fig. [3] 
we present several properties of the ground state as we 
increase the interaction, for the case of attractive interac- 
tions. The first magnitude is the population imbalance. 
As seen from the figure, it remains zero until a certain 



value (4 for N = 50) of A, then it grows abruptly and 
approaches 1 as A increases. The semiclassical approxi- 
mation would predict such behavior to occur at A = 2. 
The discrepancy between the semiclassical prediction and 
the observed quantum behavior is diminished as the num- 
ber of atoms, N , is increased. The figure also shows a 
shaded region corresponding to CTz = ^ {z^) — (z)^. In 
the region where the semiclassics fails, the dispersion of 
z becomes large. This agrees with the previous discussion 
of Figs, m and m and corresponds to the region where the 
ground state gets wider and, eventually, cat-like ^. The 
authors of Ref. [2Cl| consider the population imbalance as 
a suitable order parameter to characterize the transition. 

The one body density matrix turns out to be a good 
indicator of where the semiclassics fails to describe the 
full quantum results. In the third panel we present the 
condensed fractions, ni and n2, for different values of A. 
In the semiclassical limit, the system remains always fully 
condensed regardless of the value of A. This is clearly not 
the case in the transition region. For = 50, the exact 
dynamics differs in the region 2 < A < 4. In that region 
there are two MSPS, thus being impossible to describe 
the system within a mean field formalism. The macro- 
occupation can be traced back to the four elements of 
the one body density matrix, see panel (d). The off- 
diagonal ones are well described by the semiclassics in 
the whole considered domain, but the diagonal ones re- 
main constant (=1/2) for A < 4. The dispersion of the 
diagonal matrix elements is directly related to az , and is 
again large in the transition region. The dispersion of the 
off-diagonal elements is always of the order of 1% (thus 
explaining the agreement with the semiclassical picture) . 

The large quantum fluctuations seen in z, or pa, cannot 
be described in a mean field description. The inset in 
panel (a) of Fig. [3] shows CTz computed assuming the state 
of the system corresponds to a mean-field state ^mf = 
0))]®^, with 6,(j) taken from their semiclassical 
values. As shown, both the full quantum result and the 
mean-field one agree at A = and for A > 4, but strongly 
disagree in the transition region. Finally, the spread of 
the state in the Fock basis, S, is presented in (b). The 
function has its maximum spread for A ^ 2.4, then falls 
and has an abrupt fall off when the ground state evolves 
from cat-like to self-trapped, A '-^ 4. 

The corresponding evolution of the two MSPS is given 
in fig.m There we plot in the QL), \R)), (|+), |-)) plane 
their evolution as we increase the value of A. The exact 
results behave as follows: for < A < 2 the ground state 
of the system is fully condensed on the |-|-) state. It is 
thus symmetric, z = 0. For 2 < A < 3.5 the same two 
single particle states are macro-occupied in the entire do- 
main, see Fig. [3] -0i = |+), and V'2 = |— )• The system 



^ Employing the language of other authors, e.g. Ref. Il6ll . we would 
have that for such region the semiclassical assumption does not 
hold, ({Jz, J4> ^ 2{Jz)(Jz). 



5 




FIG. 4: Evolution of the MSPS of the ground state of the 
system as a function of A. The dashed and sohd hues depict 
the BH calculation and correspond to i/'i s.nd xp2, respectively. 
The MSPS are represented as = cos(6'/2)|L)+sin(e»/2)|ii). 
Their modulus is set to the condensed fraction of the state, 
rii. The dotted line corresponds to the semiclassical predic- 
tion. All lines are obtained by evolving from A = 0, to 15. 
BH results are given for A'^ = 25 and A'^ = 50. (This fig- 
ure corresponds to a vertical cut of the Bloch sphere used in 
Fig.H) 
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FIG. 5: Evolution of the population imbalance, thick-solid, 
and of the two macro-occupied condensed fractions, ni and 
n2, dashed and dotted lines, as a function of time for a state 
initially prepared with ^(0) = 1, *(0) = \N,0). The thin 
solid line corresponds to the semiclassical calculation of the 
population imbalance. The semiclassical values of ni and 722 
are the same for all panels and are constant and equal to 1 
and 0, respectively. The panels correspond to different values 
of NU/J, in aU cases with A^ = 50. 



in A. 



remains symmetric, z = 0. Their condensed fractions 
vary with A. For 3.5 < A < 4.5 the system has still two 
MSPS but which depart from the |±) axis. These MSPS 
change continuously as we vary A, acquiring a non-zero 
z. Finally, for A > 4.5 the system is fully condensed 
again, and ipi approaches \L) as we increase A. The fig- 
ure is constructed in such way that the departure from 
the semiclassical description becomes apparent immedi- 
ately: any point outside of the circumference of radius 
one is beyond that appoximation. The dynamics is there- 
fore mean-field like, both in the extreme Rabi regime, 
NU/J -J> 0, and in the self-trapped regime, NU/J > 4.5 
for TV = 50. The transition region cannot be described 
within a mean-field theory. The semiclassical approxi- 
mation fails to describe the observed behavior almost in 
the entire domain of the transition. 

The dependence of the described static properties as 
we vary N can be summarized as follows: increasing the 
number of atoms, the transition region gets reduced and 
thus the agreement with the semiclassical results, which 
predict a mean-field picture, is improved. As an example. 
Fig. |4] shows also the case of = 25. 

The dependence on the bias, e is qualitatively similar: 
as e is increased, the jump shifts to smaller values of A. 
In possible experimental set-ups it is however important 
to consider, as done in this work, e which are almost 
negligible to ensure the transition region is broad enough 



V. DYNAMICS: THE GROUND AND THE 
HIGHEST ENERGY STATE 

Up to now we have concentrated on describing the dif- 
ferent phenomena occurring near the transition region 
by analyzing the static properties of the spectrum of the 
Hamiltonian as we varied the parameter A. Now we will 
consider the dynamical consequences. To this extent, we 
analyze the dynamics of the system with a fixed number 
of atoms, N, with a maximal initial population imbal- 
ance, ^(0) — 1. 

The simplest quantity which shows already quantum 
fluctuations is the population imbalance. In figure [5] the 
discrepancies between the semiclassical result and the ex- 
act solutions for the time evolution of the BH model are 
easily spotted. The BH model, with finite N, always 
brings in some damping to the oscillations, followed by 
revivals, in the Josephson regime fT3|. The semiclassical 
results are seen to be accurate both for the case of no 
interaction among the atoms, Rabi oscillations, or the 
most self-trapped case. In the same figure we depict also 
the condensed fractions, ni and 712. The semiclassical 
prediction for these quantities is always ni{t) = 1 and 
"■2(0 = 0. In the finite N case however, we see how the 
condensed fractions do depart from 1 and 0, specially for 
the case 1 < NU/J < 4. 

In the case of attractive interactions the transition 
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FIG. 6: 3D representation of the eigenvectors of the one body 
density matrix, V'i(2) for NU/J — 0,1,3,7 as a function of 
time. The total time is 5tRabi. The points are represented as 
follows: 9 and (f> are obtained from, \tpi(2)) ~ cos{6 /2)\L > 
+e"^ sin(6'/2)|7?). The distance to the origin is equal to ni(2). 
The initial state is always z — 1 with complete condensation 
on the \L) state. We mark also the |±) single particle states. 

from Josephson to the self-trapped regime is completely 
related to the properties of the ground state of the sys- 
tem. Essentially, an initial state prepared with large im- 
balance and with A larger than a critical value will re- 
main trapped due to its large overlap with the ground 
state of the system. Correspondingly, the self trapping 
for the case of repulsive interactions is directly related to 
the highest excited state of the Hamiltonian, which is the 
one with larger imbalance in this case. 

In fig ini we depict the time evolution of the two eigen- 
vectors of p. The length of the eigenvector, ipi, in the 
Bloch sphere is set to its condensed fraction, n^. There- 
fore, a fully condensed evolution will remain always on 
the surface of the sphere of radius 1. All other situations, 
involving two macroscopically occupied states, would fall 
inside the sphere. 

The non-interacting case, A = 0, sets the Rabi oscilla- 
tion time, tRabi = tt/J. The population imbalance per- 
forms periodic oscillations of maximum amplitude with 
the corresponding frequency, WRabi = 2J. The system 
remains always condensed on a single particle state, see 

Fig.m 

A small interaction, NU/J = 1, already changes the 
picture. First, in this case the system is no longer con- 
densed at all times. Now, see second panel of fig [5] and 



the corresponding in fig. [6l as time increases the high- 
est condensed fraction, ni, goes from 1 to ^ 0.6 after 
t > lOtRabi. The two MSPS rotate around the (|+), |-)) 
axis, and slowly approach the |-|-) and |— ) states with a 
certain condensed fraction. 

For NU /J — A self-trapping starts. The condensed 
fraction is large, ~ 0.9, but there are sizeable fluctua- 
tions. The behavior of the MSPS is now different, in this 
case they do not approach the |-|-) axis, but instead keep 
rotating around an axis closer to the \L) vector. 

Larger interactions make the system more self-trapped 
and also to remain mostly condensed on a state progres- 
sively closer to \L). The size of the fluctuations seen in 
the condensed fractions decrease as we increase NU/ J. 

VI. CONCLUSIONS 

We have scrutinized the transition from the Joseph- 
son to the self-trapped regime in BECs, for both at- 
tractive and repulsive atom-atom interactions. First, we 
have demonstrated the impossibility of a mean-field de- 
scription of the quantum transition for a finite number 
of atoms, N, for attractive interactions. The nature of 
the transition, which involves the spontaneous breaking 
of the left-right symmetry, is governed by large quantum 
fluctuations which are not captured by a mean-field de- 
scription of the problem. The ground state of the system 
in the transition is built of two macro-occupied single 
particle states. Both for attractive and repulsive inter- 
actions we have shown how the self-trapping regime is 
related to the existence of imbalanced eigenstates in the 
spectrum. 

An extremely challenging experimental proposal em- 
anates naturally from this article, namely, the full charac- 
terization of a quantum phase transition in a BEC, either 
by considering a varying barrier height in the double- well 
or by modifying the scattering length of the atom-atom 
interaction. 
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